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ABSTRACT 


Context. The luminosity of [C ii] is used as a probe of the star formation rate in galaxies, but the correlation breaks down in some 
active galactic nuclei (AGNs). Models of the [C ii] emission from galactic nuclei do not include the influence of X-rays on the carbon 
ionization balance, which may be a factor in reducing the [C ii] luminosity. 

Aims. We aim to determine the properties of the ionized carbon and its distribution among highly ionized states in the interstellar 
gas in galactic nuclei under the influence of X-ray sources. We calculate the [C ii] luminosity in galactic nuclei under the influence of 
bright sources of soft X-rays. 

Methods. We solve the balance equation of the ionization states of carbon as a function of X-ray flux, electron, atomic hydrogen, 
and molecular hydrogen density. These are input to models of [C ii] emission from the interstellar medium (ISM) in galactic nuclei 
representing conditions in the Galactic Central Molecular Zone and a higher density AGN model. The behavior of the [C ii] luminosity 
is calculated as a function of the X-ray luminosity. We also solve the distribution of the ionization states of oxygen and nitrogen in 
highly ionized regions. 

Results. We find that the dense warm ionized medium (WIM) and dense photon dominated regions (PDRs) dominate the [Cii] 
emission when no X-rays are present. The X-rays in galactic nuclei can affect strongly the abundance in the WIM, converting 
some fraction to and higher ionization states and thus reducing its [Cii] luminosity. For an X-ray luminosity L(X-ray) > 10''^ erg 
s“' the [Cii] luminosity can be suppressed by a factor of a few, and for very strong sources, L(X-ray) >10'*'* erg s“* such as found for 
many AGNs, the [Cii] luminosity is significantly depressed. Comparison of the model with several extragalactic sources shows that 
the [C ii] to far-infrared ratio declines for L(X-ray) >10"** erg s“*, in reasonable agreement with our model. 

Conclusions. We conclude that X-rays can suppress the C^ abundance and, therefore, the [C ii] luminosity of the ISM in active galactic 
nuclei with a large X-ray flux. The X-ray flux can arise from a central massive accreting black hole and/or from many smaller discrete 
sources distributed throughout the nuclei. We also find that the lower ionization states of nitrogen and oxygen are also suppressed at 
high X-ray fluxes in warm ionized gas. 
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1. Introduction 


The fine structure line of ionized carbon, [Cii], at 158 pm 


(IStacev et alJl2010 

; Ide Looze et alJl201 ll iDe Looze et alJl2014 ; 

Pineda et alJ 12014 

), as is the far-infrared dust emission (e.g. 

Malhotra et al. 

|1997 

iLuhman et alJ 119981: IStacev et al.1 1201 Ol 


De Looze et al.ll2()14 ). In principle the energy output of bright 


young stars is transferred to the gas and dust and reemitted ei¬ 
ther in the important far-infrared cooling lines such as [Cii] 
or as continuum dust emission. In extragalactic sources the lu¬ 
minosity of [C ii] is often correlated with the far-infrared dust 
emission thus reinforcing the suggestion that these trace the 
SFR. However there are some notable exceptions in which 
the luminosity of [Cii] is significantly suppressed with re¬ 
spect to the infrared luminosity. The relationship between 
the intensity of [Cii] and other sta r formation tracer s ap¬ 
pears to work in the Gala ctic disk (iPin eda et all 1201 4t) and 
most extragalactic sources (iDe Looze et al. 2014 ). however it 
breaks d own somewhat in many luminous infrared galaxies 
(LIRGs) (iDfaz-Santos et al. 20141) and ul t ra luminous infrared 
galaxies (ULIRGs) (Malhotra et al.l 119971 : iLuhman et al.lIl998t 
iRigopoulou et alJ 1201 4l) an d even more so in many active 
galac t ic nuclei (AGN) (c.f . IStacev et alJ 120101: ISargsvan et alJ 
l2014t iGullberg et alJlSoi^ . The assumption is that there is a 


deficit in [C ii] lumi nosity in these sources and in particular in 
their central regions (Malhotra et al.l 1997 : Luhman et al]|1998l : 


iRigopoulou et al.ll201^Wilbers et alfcOlSl).' 


Several explanations have been offe red for the decrease in 
the [Cii] to fardnfrared luminosity (e.g. Nakagawa et akl 


Malhotra et alJ 1199'^ iLuhman et akl Il998r Stacey et al 


1995 


201C: 


Gullberg et al. [201^ . including: 1) the [Cii] emission is opti 


cally thick and saturated, if it arises from very dense photon 
dominated regions (PDRs), where the excitation temperature of 
the ^P 3 / 2 -^P 1/2 transition is thermalized and does not radiate as 
efficiently as the dust with increasing energy input; 2) for high 
ratios of UV flux to density the dust grains are positively charged 
and are less efficient at heating the gas vi a photoeiected elec¬ 
trons, thus the ^P^n level is less populated (IWolfire et al.lll990l : 
iKaufman et alJ[7998h : and, 3) a reduction in the column density 
of ionized carbon in the PDRs due to a soft ultraviolet radiation 
field that can heat the dust but not maintain a large layer of C*^. 
However, detailed models of the source of [C ii] are lacking and 
none of the proposed mechanisms appear to explain the reduc¬ 
tion in all sources. 


In this paper we propose that X-rays can also produce a 
deficiency in the [C ii] to far-infrared (FIR) luminosity ratio in 
the presence of a large X-ray flux at energies above ~1 keV. 
Under these conditions, a significant fraction of carbon is in 
higher ionization states (C^*^, etc.) in some of the regions 
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responsible for [C ii] emission, in particular the highly ionized 
ISM gas components. We also propose that X-ray photoion¬ 
ization from strong X-ray sources can reduce the abundance 
of the lower ionization states of nitrogen and oxygen in the 
ionized gas and therefore the emission from their far-infrared 
hne-structure lines. For example, the [N iii] 57 fim line is ob¬ 
served to be comparable to the [N ii1 1 22 /rm line in extragalac- 
tic sources (iGracia-Carpio et Sl201 Ih . The abundance of 
will be comparable to or greater than that for because the 
critical density for exciting the hne-structure leve ls is much 
higher than that of those for N^ dMadden et al.ll2013h . Previous 
work on the effects of X-rays on chemistry have focused on the 
X-ray Dissoci a tion Regions (XDRs) of molecular clouds (e. 


Langer 19781 Krolik & KallmanI 119831: Maloney et al. 199t 
Meiierink & Spaansll2005l) or YSOs ( Bruderer et al. 20091) and 


protoplanetary disks (lAdamkovics et al]l201 llf*l 

Here we consider the effects of high X-ray fluxes on the 
interstellar gas environments that emit in [Cii]. These include 
primarily the highly ionized gas, which is composed of sev¬ 
eral components including the warm ionized medium (WIM), 
Hot Ionized Medium (HIM), the dense ionized skins surround¬ 
ing clouds (DIS), and Hii regions. A large X-ray flux will also 
affect the [C ii] emission from diffuse atomic hydrogen clouds 
and, to a much lesser extent, the photon dominated regions 
(PDRs) surrounding molecular clouds. We also consider the ion¬ 
ization balance of nitrogen and oxygen in the WIM. The em- 

S ihasis in this paper is different from that in the XDR models 
Maloney et al.l[l9^ iMeiierink & Spaansll2005l) in that we do 
not evaluate the molecular abundances in the dense UV shielded 
molecular clouds, but focus instead on the abundance of the car¬ 
bon ionization state in the major ISM components. Whether X- 
rays play a signihcant role in reducing the abundance of (and 
other gas tracers such as N^, 0° and in the ISM and, 

thus, [C ii] emission will depend on details of the distribution of 
X-ray sources, their luminosity, and the properties of the ISM. 
In addition, a high flux of X-rays can compete with UV he ating 
result ing in an increase in the far-infrared (FIR) luminosity (IVoid 
Il99lh and therefore further reducing the [C ii] to FIR ratio. 

We will show that the X-ray flux can reduce the abun¬ 
dance of primarily in highly ionized regions and that, un¬ 
der the right conditions, will reduce the [C ii] emission signif¬ 
icantly in galactic nuclei. In the Galactic disk th e WIM con- 
tributes about 20% to 30% of the [C ii] luminosity (iPineda et al.l 
1201 3t IVelusamv & Langerll2014 ) using the scale heights der ived 
from rC III dUanger et al.ll2014a : IVelusamv & Langeill2014h . In 
the Galactic Central Molecular Zone (CMZ) the WIM densities 
are much higher and may contribute an even larger percentage 
of the [C ii] luminosity. We conclude that this mechanism is im¬ 
portant in some active galactic environments with large X-ray 
luminosity and must be included along with the previously sug¬ 
gested mechanisms to explain the reduction in the [C ii] to far-lR 
luminosity in galaxies. 

X-rays are an important radiation component of active galac¬ 
tic nuclei (AGNs) composi n g up to 40% of of th e total luminos¬ 
ity dMushotzkv et al.l[l993h . Ebrero et~^ (l2009l) summarize the 
soft (0.5 - 2 keV) X-ray luminosities in a large sample (> 1000) 
of AGNs (reproduced in Fig. [TJ whi ch have lumino s ities i n the 




wasawa 


etal.N2011h in ; 


range 10'^*’ erg s“* to 2x10"^^ erg s^ 
recent study of Luminous Infrared Galaxies (LIRGs) found that 
43 out of 44 sources had X-rays luminosities ranging from ~10^*^ 


* There is also a body of work on the carbon ionization balance in 
very hot di ffuse ionized gas due to electron collisional ionization of 
carbon fe.g. lGnat & SternberdUOPTh . 



z 


Fig. 1. Soft X-ray (0.5 - 2.0 keV) luminosity for a sample of 
1009 AGNs a s a fun ction of red shift z reproduced from Fig. 1 in 
lEbrero et al.l (l2009h . The inse rt key identif i es the X-ray survey 
used to compile the data (see lEbrero et al.l (l2009h for their def¬ 
initions). The soft X-ray luminosities for this sample of AGNs 
range from ~10"^‘’ to ~2xl0"^’ erg s“'. 


erg s“* to ~10"^^ erg s“* in the 0.5 to 10 keV range. The mean 
radius of the soft X-ray emission in their sample was 5.3 kpc 
and so represents emission throughout the LIRG, however they 
found that the distribution of the hard X-ray emission is much 
more compact. X-ray luminosities greater than 10"^° erg s“' rep¬ 
resent an important ionization source that must be taken into ac¬ 
count in understanding the conditions of the interstellar medium 
and the corresponding emission in far-infrared lines from ions. 

Even in more quiescent galaxies X-rays are important in the 
nucleus. Eor example, the Galactic CMZ is observed to have ex¬ 
tended d iffuse X-ray emiss ion as well as numerous X-ray point 
sources (Muno,et_ay200^, including an accreting stellar black 
hole (iHeindl et al.l Il993lk In addition, the massive black hole 
near SgrA is a potential X-ray source. It has been suggested that 
X-rays could provide the > 10 s ' ioniza tion rate of hydro¬ 
gen needed to explain the H 3 ^ observations of iGoto et alJ(l2014l) . 
The cross section for ionization of heavy atoms by X-rays is or¬ 
ders of magnitude larger than hydrogen, primarily due to K shell 
ionization at X-ray energies greater than a few hundred eV. In 
addition. X-ray photoionization of hydrogen produces energetic 
photoei ected electrons which produce secondary ioniza tion of 
carbon (iMalonev et alJ[T9^ IMeiierink & Spaan^l2005l) . Thus, 
X-rays can play an important role in ionizing carbon and other 
metals. The penetration depth of 1 keV X-rays in the diffuse ion¬ 
ized gas in the CMZ is a column density V(H^) > fewxlO^^ 
cm“^ and in the neutral gas, is V(H-i- 2 H 2 ) ~ several xlO^* cm“^. 
Thus X-rays suffuse galactic nuclei and X-ray photoionization 
needs to be considered in the ionization balance of metals in the 
interstellar medium. 

Our paper is organized as follows. In Section |2] we develop 
the reactions and ionization balance model, while in Section [3] 
we present sample model calculations of the distribution of ion¬ 
ization states of C, N, and O. In Section |4] we model the [C ii] 
luminosity as a function of X-ray luminosity, and discuss the 
consequences and implications of X-ray photoionization on the 
[Cii] emission from active galactic nuclei. In Section|4]we also 
validate our [C ii] emissivity model for galactic nuclei by calcu- 
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lating the luminosity from the CMZ and comparing it to obser- where 
vations. Finally, we summarize our results in Section|5] 


2. X-ray ionization modei for carbon 

In this section we develop the X-ray photoionization balance 
equations for the ionization states of carbon. Multiply ionized 
carbon can be produced by L- and K-shell photoionization by 
X-rays. K-shell photoionization dominates L-shell photoioniza¬ 
tion at energies above its threshold, Ek ~ 0.3 - 0.4 keV for car¬ 
bon. Furthermore, these higher energy X-rays penetrate further 
in the ISM than the low energy X-rays that dominate L-shell 
photoionization because the L-shell photoionization cross sec¬ 
tions are large for hydrogen and carbon near their thresholds ~ 
11 to 15 eV. We consider only X-rays with energies > 1 keV 
because they have large enough cross sections to influence the 
ionization balance but not so large that they cannot penetrate 
large column densities. 


2.1. Carbon ion balance equations 


Here we describe the reduced set of reactions that dominate the 
ionization balance of carbon in the presence of X-rays. The X- 
ray photoionization processes are: 


hvx + -H e 

hvx + -H 2e 


( 1 ) 

( 2 ) 


where j labels the initial ionization state, CL and ranges from 
0 to 5 for L-shell ionization (Equation [TJ and 0 to 4 for Auger 
K-shell ionization (Equation |2]i. The X-ray photoionization rate 
from is designated and for L-shell and K-shell, re¬ 
spectively. 

The electron recombination reactions are, 

C<Ai)+ + g ^ H- hv, (3) 


where the reaction rate coefficient from j + 1 to j is labeled ki 
Ions can charge exchange with atomic hydrogen, 

+ H ^ . (4) 


where the reaction rate coefficient from j + 1 to j is labeled kl^. 
In regions with molecular hydrogen generally reacts rapidly 
via ion-molecule reactions. 


C0+1)+ + H 2 


(5) 

C0+1)+ -1- H 2 

+H* +H^ 

(6) 

C0+1)+ H- H 2 

-H FU -H // 

(7) 


except for which reacts slowly with H 2 via radiative associa¬ 
tion, 

+H 2 ^ + hv, (8) 

we label these ion molecule reactions . To calculate the num- 

im 

her density of carbon ions we need to solve a set of balance equa¬ 
tions, 

^XK ) “ Lx/r 


^^■ = k>(e)+ki.xn(H) + kJ^n(H2). 


(10) 


for i -Oto jmax-^ (note that ^1'"" = 0 and 0 in Equation^ 
because there is no ionization out of and no recombination into 
the maximum charge state, respectively, and that = 0 be¬ 
cause there is no recombination from the lowest charge state). In 
addition we need the number density conservation equation. 


J]nj(Cn^n,(C) 


0 


( 11 ) 


where n,(C) is the total number density of all carbon charge 
states. Eor the galactic center interstellar medium the X-ray en¬ 
ergy range of interest is greater than 1 keV and we can neglect 
the L-shell photoionization rate (see Section ITTI) . These equa¬ 
tions can be solved analytically and used to calculate the frac¬ 
tional ionization of each species, fj = rijln,, where i labels the 

ionization state, and 2 /i' - 1 ■ 

0 

2.2. X-ray photoionization cross sections 

IVerner et al.l (Il993h provide fits for X-ray photoionization cross 
sections for C° to for the Is-shell (K-shell), C** to for 
the 2s-shell (L-shell), and C** and for the 2-p shell (L-shell). 
The 2s- and 2p-shell cross sections contribute less than 5% to the 
total cross section above 0.4 keV. To facilitate solving for we 
fit tfie cross sectio ns for X-ray photoio nizat ion of C°, C^, 

a nd from IVerner et al.l (Il993h and IVerner & Yakovle^ 
(dial with a power law fit above 0.4 keV. The fits are very sim¬ 
ilar for all the K-shell reactions and have a form, 

cr(E) - croE^^y, and for simplicity we adopt an average value, 
cTo = 4.1 X 10“^° cm^ and b - 2.8. In SectionOwe will combine 
cr(E) with the X-ray photon flux to calculate the X-ray photoion¬ 
ization rates. As we are only interested in solving for the effects 
of X-rays on the carbon ionization balance in regions where neu¬ 
tral carbon is readily ionized by UV we will neglect the balance 
equations for C°. This assumption is valid in the presence of UV 
which ionizes neutral carbon rapidly in unshielded regions in the 
interstellar medium, otherwise there would not be any and it 
would recombine to C°. Eurthermore, the ionization of C° by 
UV is generally orders of magnitude larger than that by X-rays. 


2.3. Electron recombination reaction rate coefficients 

The reaction rate coefficients for radiative and dielectronic re¬ 
combination of electrons with carbon ions up to have been 
calc ulated bylNahar & Pradh^ (Il997h and iBadnell et al.l (l2003h 
and iBadnelll (l2006l) with radiative recombination dominating 
at low temperatures (as in diffuse clouds) but dielectronic re¬ 
combination becoming more important at the higher tempera¬ 
tures found in the WIM and DIS. We use the values for the 
total recombination rate coefficients published in Table 5 of 
lNahar&Pradl^(ll997h in our model calculations. The reac¬ 
tions are given in Table[T]and we list the reaction rate coefficients 
for two characteristic gas temperatures, 100 K and 8000 K. 
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Table 1. Electron recombination reactions Table 3. Reaction rate coefficients of C"^ with H 2 


Reaction" 

lOOK 

8000K 

Reaction 

kilOOKf 

C+ -H e ^ C" + hv 

9.1xl0-‘" 

7.3x10-'" 

C+ -t- H2 ^ CHJ + hv 

6x10-'° 

+ + hv 

4.2x10-“ 

6.9x10-'^ 

-H H2 ^ C° -H H+ H+ 

1.2x10-'° 

-1- e ^ + hv 

3.2xl0-‘° 

1.8x10-'' 

+ H2 ^ C+ + H+ + H+ 

1.2x10-'° 

-1- e ^ C 3 + + hv 

1.6x10-'° 

1.0x10-“ 

C‘'+ + H2 ^ C 2 + + H+ 4 - H+ 

1.2x10-'° 


a. The reaction rate co efficients, in cm^’ s *, are from Table 5 of a. In units of cm^ s“'. 
Nahar & Pradhanl j 19971) . 


2.4. Charge exchange with atomic hydrogen 

Reaction rate coefficients for charge transfer between carbon 
io ns and atomic hydrogen hav e been summarized for to 
bv iKingdon & Fertandl (1 199611 afong with fits in the form, 

kcx{T 4 ) = aT^^[l + ce‘^T4] cm^s-' (12) 

where T 4 is the kinetic temperature in units of 10"^ K, and the 
fitting coefficients are given in their Tabfe 1. The charge ex¬ 
change of with H is very smaff owing to the iarge endother¬ 
mic barrier, AE~ 2.5 eV (due to the different ionization poten- 
tiafs of C and H). Typicaffy most exothermic charge exchange 
reactions of highfy ionized metafs proceed rapidfy under ISM 
conditions with a reaction rate coefficient of order 10“^ cm^ s *, 
but that is not the case for doubfy ionized carbon. Afthough 
charge exchange of with H is energeticaffy favorabfe, its 
cross section is smaff owing to the nature of the crossings in the 
intermediate motecutar p otentiaf dMcCarrott & Vafironl Il97-5t 
iKingdon & Eertandl 19961) . However, the higher ionization states 
of carbon do charge exchange rapidfy. In Table |2] we list the re¬ 
action rate coefficients for carbon ions at 8000 K and 100 K, 
typical of temperatures found in the warm ionized medium and 
diffuse atomic hydrogen gas, respectively. 


Table 2. Charge exchange with H reaction rate coefficients 


Reaction 

fcrx(8000 K)" 

fecx(100 K)" 

C+ + H C" + H+ 

-7x10-'" 

«io-"° 

C2+ 4- H ^ C+ H+ 

l.lxl0-'2 

1.3x10-'° 

C"+ 4- H ^ C2+ 4 - H+ 

3.1x10-° 

1.5x10-° 

4- H ^ C"+ -H H+ 

2.1x10-° 

2.8x10-° 


a. In units of cm^ s *. 


2.5. ion moiecuie reactions with H 2 

reacts slowly with molecular hy drogen via associativ e re¬ 
combination, -I- H2 —» CHj + hv (iMcElrov et al.ll2013h . and 
this reaction does not play a significant role in the carbon ion¬ 
ization balance outside of dense PDRs. In contrast, the reactions 
of C"^ with H 2 for n > 2 are likely to be fast, ~10 to 10“® 
cm^ s“' for most cases (lLangerlll978h . Thus the admixture of 
a small amount of molecular hydrogen can haye a large effect 
on the carbon X-ray photoionization balance. The only multiply 
charged carbon ion molecule reaction that h as been measured in 
the lab is, + H 2 (ISmith & Adamsll98lh and it has a reaction 
rate coefficient of ~ 1.2x10“'° cm^ s“* with the dominant chan¬ 
nel ^ -H -H C comprising ~90% of the total. We assume 
that higher charge states haye the same reaction rate coefficient 
and end products, although the actual product channel is not crit¬ 
ical for the model, and these are listed in Table [3] 


3. Results 

The gas in galactic nuclei is exposed to a different radiation en- 
yironment than that in the galactic disks with intense UV and 
X-ray radiation fields due to the presence of young O-type stars, 
supernoyae, stellar black holes, and a massiye black hole at the 
center. In AGNs this enyironment reaches an extreme with X- 
ray luminosities, L(X -ray), in excess of lO""* erg s“' and up to 
~2xl0"'^ erg s“* (e.g. lEbrero et al.ir2009l) . These large luminosi¬ 
ties are sufficient to influence the distribution of higher ionized 
states of carbon in galactic nuclei and depends on the flux of 
soft X-rays with energies aboye the K-shell ionization threshold 
(~0.3 - 0.4 keV). If the source of X-rays in the galactic center is 
due to an accreting black hole the X-ray luminosity flux. 




L(X-ray) 


(13) 


In the notation of iMaloney et al.l (Il996h . Fx ~ 8.4xl0%4rp,^ 
erg cm“^ s“', where L 44 - Lx/10^ erg s“* and rpc is the distance 
to the X-ray source in pc. Thus, neglecting absorption aboye 1 
keV, Fx can be yery large throughout most of a galactic nucleus 
and these X-rays will haye a major impact on the properties of 
the neutral and ionized gas. While X-ray emission from massiye 
black holes can effect the ionization balance, chemistry, and ther¬ 
mal properties oyer large yolumes of the galactic center, other 
weaker sources can haye comparable influence on nearby clouds. 
For example, strong supernoya shocks can produce luminosities 
Lx ~ 10^° to 10^° e rg s ' and stellar b lack holes can produce 


-10^’ - 10^ ^ erg s ' (iHeindl et al.l[T9^ . 


iMaloney et al.l d 19961) modeled the effects of enhanced UV 
and X-ray fluxes on the thermal and chemical balance of molec¬ 
ular gas in galactic nuclei, generating a grid of XDR models that 
coyered a range of X-ray radiation field energy fluxes Ex = 0.1 
to 10° erg cm“^ s“'. They set a lower limit of 1 keV to guarantee 
that the X-rays will penetrate sufficient depth in to the molecular 
cloud to affect the chemical and thermal balance. [Meiierink et al.l 
(I2OO7I) also modeled the effects of enhanced UV and X-rays on 
molecular gas in galactic nuclei, but coyered a narrower range of 
energy flux, Fx = 1.6x10“^ to 1.6x10^ erg cm“^ s“'. Here we 
consider the effects of large X-ray fluxes on the carbon balance 
in the ionized ISM, the diffuse atomic hydrogen clouds, and the 
photodissociation regions (PDRs) of dense molecular clouds in 
galactic nuclei. We begin by calculating the X-ray photoioniza¬ 
tion rates. 


3.1. X-ray photoionization rates 

Here we calculate the X-ray ionization rates for carbon ions pa¬ 
rameterized in terms of the X-ray energy flux, Fx. The X-ray 
photoionization rate for carbon ions of charge state j is 

; r^"'“ dJ 
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where dJIdE is the photon spectral distribution of X-rays as a 
function of energy in units of photons cm“^ s ' keV *, and Emin 
and Emax are the minimum and maximum X-ray energies, re¬ 
spectively. The spectral distribution is usually expressed as a 
power law distribution, dJ/dE = For Seyfert galax¬ 

ies r is typically in the r ange 1.4 to 2.0 below 1 00 keV with a 
distribution average ~\.l dMushotzkv et al.lll993h. However, the 
spectr al index can vary with energy, for example [iwasawa et alJ 
( 1201 Ih found different spectral slopes for hard a nd soft X-rays 
in a s ample of LIRGs. For active ga lactic nuclei dM a onev et^ 
Il996ll adopted F ~ 2 for AGNs and iMeiierink et akT i 2007 ) F = 
1.9 for the X-ray spectrum of an accreting massive black hole. 
Here we will adopt F = 1.9. We will set Emin - 1 keV because, 
as discussed above, lower energy X-rays can be neglected except 
very near the X-ray source, and Emax is usually taken to be 100 
keV. 

To calculate the photoionization rate as a function of the en¬ 
ergy flux we will vary 7o to represent different X-ray fluxes. We 
rewrite the spectral energy distribution as. 


dJ 


FxUl) 

pr 

^keV 


(15) 


where Jq(1) is the value that produces a luminosity flux of 1 erg 
cm“^ s“' for a given spectral index F. We And 7o(l) by solving 




£ 


E--—dE = 1 (erg cm s 

pT ° 

^keV 




(16) 


substituting F =1.9, integrating from 1 to 100 keV, and con¬ 
verting units from keV to ergs. This conversion yields, 7o(l) 
= 1.1x10^ photons cm“^ s keV^'. We now express dJ/dE 
- 1.1x10''F'x/£’''^ where Ex is in units of erg cm“^ s“*. 
Substituting dJ/dE and cr{Ekev) into EauationfT4lvields 

r""' dE 

= 1.1 X 10'F;frro (17) 


This equation is valid at the boundary of the cloud as long as 
we can neglect absorption of the X-rays. For 1 keV, our assumed 
lower threshold, where most of the ionization occurs, the absorp¬ 
tion opacity is less than one for a hydrogen column density < (5 
- 10)xl0^' cm“^, depending on the metallicity. Only the GMCs 
have sufficient column density to absorb X-rays at this energy, 
however, as discussed the X-ray ionization has little effect on 
the [C It] emission from the molecular hydrogen boundaries of 
the cloud and therefore we neglect the attenuation of X-rays in 
determining the luminosity from this layer. Substituting ctq, F, 
and b into Equation[T7]and integrating from 1 to 100 keV yields, 
1.2x10-'2Fx s-'. 

In the neutral gas clouds the electrons ejected by X-ray pho¬ 
toionization primarily of hydrogen produce secondary ioniza¬ 
tion. The multiplicative effects of these secondary ionizations 
of metals is of order a few and here we adopt a factor of tw o 
for carbon (see lMalonev et aO 19961 : IMeiierink & Spaansll2005h . 
However, in the ionized WIM and DIS the photoejected elec¬ 
trons scatter via Coulomb interactions off the ambient electrons. 
This scattering is large co mpared with the el ectron ionization 
cross sections for C"^ (cf. ISuno & Katoll2003) and the photoe¬ 
jected electrons are rapidly thermalized. Therefore above a frac- 
tional ionization of a c ouple of percent, thermalization is so rapid 
dMalonev et alJlT99^ that the probability of secondary ioniza¬ 
tion of carbon can be neglected in the WIM and DIS. 


The X-ray photoionization rate of carbon ions can be quite 
large in galactic nuclei and this process will result in the conver- 
sion of to higher io nization states in diffuse ionized regions. 
IMeiierink et al.N2007h considered Ex = 1.6x10“^ to 1.6x10^ 
erg cm~ ^ s ' in their XDR m odels of dense gas in galaxy nu¬ 
clei and iMalonev et alJ (1 19961) studied a range of 0.1 to 10® erg 
cm“^ s '. In this Section, for illustrative purposes, we will usu¬ 
ally consider an X-ray flux of Ex - 10“^ to 10"' erg cm“^ s“', 
but sometimes consider higher values. The lower limit reflects 
an ionization rate at which X-rays are unimportant in the ion 
balance. In the following subsections we present the distribution 
of ionized states of carbon as a function of X-ray flux, under 
different conditions of electron, atomic hydrogen, and molecu¬ 
lar hydrogen abundances, corresponding to different AGN ISM 
environments. To simplify the solution we include an ionization 
state only up to C"'^, as our main purpose is to study the conver¬ 
sion of to higher states. Thus the fractional abundance /(C"''^) 
represents all charge states z >4. 


3.2. Highly ionized regions 


There are several environments in which hydrogen is completely 
ionized. One is the warm ionized medium, where the densities 
in the disk are usually low, n(e) ~ fewxlO^^ cm“^ (iHaffner et al.l 
l2009l) . but are much higher in galactic nuclei. In the Centra l 
Molecular Zone, ~ 10 cm“^ dCordes & Lazioll2003l : lRovll2013l) . 
but may be higher in AGNs. Other highly ionized regions are 
those associated with Hit sources and Hii-like dense ionized 
skins (DIS) around dense molecular clouds, where the electron 
densities may be h igh, n(e) ~ 1 - 100 cm“^ (lOberst et al.ll201 ll 
iLanger et akl 120151) . Einally, there is the hot ionized medium 
(HIM) which is detected via thermal X-ray emission and where 
hydrogen and helium are completely ionized. The HIM would 
have Tkin ~ 10® - few xlO® K and n(e) < fewxlO^® cm ® (ICoxI 
12005^ l^rriere et ^l2007l) . In general the densities for the HIM 
are too small for this component to contribute much to the [C ii] 
emission from galactic nuclei and we neglect it here. 

In Eig. |2la) we show the solutions for a low density com¬ 
pletely ionized gas, n(e) = 0.01 cm~® (fully ionized, n(H) = 
«(H 2 ) = 0 cm ®) and Tkin = 8000 K. It can be seen that as the 
X-ray energy flux increases the is first converted to then 
and at Ex ~10 ' erg cm"^ s ' becomes negligible and all 
the carbon is in the higher ionization states. For Ex > I erg cm“^ 
s ' carbon is essentially in the highest charge state considered 
here, C"'^. For galactic centers and Hii-like regions where the 
electron densities are higher, we show results for n(e) = 1 cm^® 
and 10 cm ® in Figs.|2jb) and (c), respectively. Increasing the 
electron density shifts the transition from to higher ionization 
states to larger Ex as a larger X-ray radiation flux is needed to 
offset the more rapid electron recombination to lower ionization 
states. For all the electron densities considered here there is a 
threshold in Ex where goes to zero in the highly ionized gas. 
Adding a small amount of atomic hydrogen, n(H)~ 0.01 cm 
does not have much of an effect on the distribution of ionization 
states of carbon (not shown). For the hot ionized medium the 
combination of a lower electron recombination rate coefficient 
at the higher temperatures and the lower electron density pushes 
the transition from to to low values of Ex', for example, 
for Tkin = 10® K and n(e) = 0.003 cm“® we And /(C"'") <0.1 for 
the lowest flux considered here Ex = 10 ® erg cm“^ s '. 
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Fig. 2. Fractional ionization state of to versus X-ray lu¬ 
minosity flux, Fx, for n(e)= 0.01 cm“^ (top), 1.0 (middle), 
and 10.0 cm^^ (bottom) at rkm = 8000 K, assuming a fully ion¬ 
ized region (n(H) = «(H 2 ) = 0 cm“^). 


3.3. Atomic hydrogen clouds 

We next address the question whether a large flux of X-rays will 
modify the distribution of in low density diffuse atomic hy¬ 
drogen clouds. The presence of even modest amounts of atomic 
hydrogen will keep carbon from reaching highly ionized states 
because of the efficiency of charge exchange, except for 
which does not charge exchange efficiently. In Fig. [3a) we show 
the fraction of and higher ionization states for conditions in 
a diffuse atomic hydrogen cloud with Tan =100 K, a hydrogen 
density n(FI) = 25 cm“^, a small admixture of electrons, and no 
molecular hydrogen. As can be seen in this figure only the two 
lowest carbon ionization states are present, and the higher charge 
states have very small fractional abundances because they charge 
exchange rapidly with H to lower ionization states. For Fx < 10 
erg cm“^ s“^ carbon is mainly singly ionized, however in the 
presence of an X-ray energy flux > 10 erg cm“^ s ' X-ray K- 
shell photoionization converts rapidly into which then 
charge exchanges to Thus under high X-ray flux conditions 
low density atomic hydrogen clouds would have a low fraction 
of but a high fraction of and [C ii] emission would be re¬ 
duced. In Fig.|3b) we show the same calculation for a hydrogen 
density n(FI) =100 cm^^ and we see it has the same behavior as 
the lower n(H) case but the transition occurs at a higher flux Fx. 

The results in Fig. [3] assume a very small, but non-neglible, 
electron abundance because it will be present in atomic clouds 
due to ionization of hydrogen by cosmic rays and X-rays. Under 
typical ISM conditions found in the Galactic disk the elec¬ 
tron fraction in diffuse atomic clouds is small («(e)/n(FI)<10“^). 
However, the electron fraction increases with increasing X-ray 
flux. We can estimate the electron fraction by balancing X- 
ray ionization with electron recombination. We And a fractional 
electron abundance ~10“^ for Fx - 10^ erg cm“^ s ' in a cloud 
with n(H) =10^ cm“^. However, the addition of a small frac¬ 
tional abundance of electrons (< 0.1) does not change the distri¬ 
bution of in the diffuse clouds. Finally, we note that for very 
large X-ray energy fluxes, Fx > 10^ erg cm“^ s ' and for n(H) < 
100 cm“^ that the fractional abundance of H^, /(H^) >0.1 and 
the assumption of a neutral atomic hydrogen cloud is no longer 
valid. 


3.4. Molecular hydrogen clouds 


Highly ionized states of carbon react rapidly with molecular hy¬ 
drogen and so the main reserv oirs of carbon in molecular c loud 
PDR and XDR regitm s is C^. iMeiierink & Snaanll (l200-5h and 
iMeiierink et al.l (l2007h included the produc tion of by X-ray 
photoi onization in their XDR models (see IMeiierink & Snaani 
(l200-5h Appendix D.3.1). However, as discussed above, this rapid 
reaction reduces the abundance to very low values in their 
models, even at their maximum energy flux, Fx - 160 erg cm“^ 
s“'. We illustrate the influence of a small admixture of molecular 
hydrogen, n(H 2 ) = 10 cm“^, or a 10% admixture, on the carbon 
ionization balance in a diffuse hydrogen cloud with n(H) =100 
cm“^ and Tkin = lOOK. As seen in Fig. |3c) the crossover from 
to occurs at a higher value of Fx than in a cloud with 
only atomic hydrogen (Fig.|3b)). 

The rapid charge transfer with molecular hydrogen has a 
more profound effect in dense PDRs where the H 2 suppresses 
significantly the population of higher carbon ionization states as 
is illustrated in Fig.|4|for n(H 2 ) = 10^ cm“^, where the hydrogen 
density is held constant (which may not be vali d at high X-ray 


fluxes as disc ussed below). In the XDR models o fiMalonev et al 


(1 19961 1 19971) and IMeiierink & Snaani d2005l) : IMeiierink et al 
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Fig. 3. Fractional ionization state of through versus X- 
ray luminosity flux, Fx, for a diffuse atomic hydrogen cloud with 
n(H) = 25 cm“^ (panel a) and 100 cm~^ (panel b), assuming n(e) 
= 0.1 cm~^, and rkm = 100 K and no H 2 , panel (c) includes a 
10 % admixture of H 2 . 


(I 2 OO 7 I 1 the molecular hydrogen density n(H 2 ) >10^ cm“^ and 
molecular hydrogen is efficient at suppressing the highly ion¬ 
ized states of carbon to such an extent that is necessary to go to 
very large X-ray fluxes, Fx >10^ erg cm“^ s“' to see any sig¬ 
nificant decrease in C^. The results are relatively insensitive to 
the choice of kinetic temperature because the charge exchange 
reactions with H 2 dominate the ionization balance in the PDRs . 

In Fig.|4]the hydrogen density was held constant to illustrate 
the effects of the reactions of the multiply ionized carbon with 
the ambient gas. However, the X-rays not only ionize carbon 
but also ionize the hydrogen molecules, producing i/J which 
can subsequently be destroyed in reactions with H 2 leading to 
atomic hydrogen (HJ -F H 2 —» H 3 - 1 - H, followed by reactions 
of with atoms, trace molecul es and electrons that le ad to ad¬ 
ditional production of H atoms'). iMalonev et al.l dl996h estimate 
that every direct ionization of H 2 leads to three H atoms. The 
formation of H 2 takes place on grain surfaces and the forma¬ 
tion rate depends on the grain cross section, sticking coefficients, 
and recombination efficienc y as summarized in Appendix D. 1 of 
iMeiierink & Snaani (l2005h . Typical formation rates are of order 
(1 - 3)xl0“'^n(H)n,o,, where is the total number of hydro¬ 
gen nucleons (= n(H) -F 2 n(H 2 )). Thus the PDR layers of dense 
molecular clouds become atomic hydrogen (and eventually H^) 
above a threshold Fx- We estimate that for dense PDR layers 
the fraction of H, /(H) is of order one for Fx > 10^ erg cm“^ 
s / Detailed models of PDR/XDR layers show this transition 
(see IMalonev et al]ll996t IMeiierink & SDaan^l2005h . Thus, al¬ 
though X-rays can, in principle, produce multiply ionized carbon 
ions in the H 2 gas, under realistic conditions in dense molecular 
clouds this situation is unlikely to be maintained as the molecu¬ 
lar hydrogen is dissociated at a lower energy flux than needed to 
sustain highly ionized states of carbon. As discussed below this 
process results in the PDR/XDR layers occurring deeper in the 
cloud where the X-ray flux has been attenuated. 

3.5. X-ray ionization of nitrogen and oxygen 

In addition to carbon, atomic oxygen, and nitrogen and oxygen 
ions have fine-structure transitions that are important diagnos¬ 
tics of the ISM in galaxies and are used as tracers of the star 
formation rate (see, for example IPe Looze et al.l 1^141) . These 
atoms and ions also play an important role in the thermal bal¬ 
ance of the galactic gas. The O and N atomic and ionic lines 
also show intensity deficits (e.g. iLuhman et al.ll 199^ iDale et all 

l2004t l&acia-Carpio et al.ll201 Ih. 

In this paper we do not model the complete ionization distri¬ 
bution of these species and their far-infrared emission but instead 
focus on a few key differences in their chemistry in highly ion¬ 
ized gas as it is likely that much of the emission from galaxies 
come from these ISM components. These differences alter the 
distribution of their higher ionization states as compared with 
carboifl First, the ionization potential of N (14.534 e'V) is sig¬ 
nificantly greater than hydrogen (13.598 e'V) so photoionization 
occurs only for EUV photons somewhat beyond the Lyman limit 
(911.8 Angstrom). In contrast the ionization potential of oxygen, 
13.6181 eV, is only slightly above that of hydrogen but still be¬ 
yond the Lyman limit. In the warm ionized medium charge ex¬ 
change reactions with protons are able to sustain a significant 
fraction of and even in the absence of LUV or X-rays. In 
the case of oxygen the energy difference is small (in temperature 

^ An example of the differences in the distribution of C, N, and O 
ionization states in the presence of X-rays can be seen in the models of 
lAdamkovics et al.l ll201lh for protoplanetary disks. 
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Fig. 4. Fraction of carbon ionization states versus X-ray lumi¬ 
nosity flux, Fx, for a dense molecular hydrogen gas, n(H 2 ) = 
10^ cm“^ and Fkin =100 K, with an admixture n(H) = 1 cm“^ 
and n(e) = 0.1 cm“^. Carbon remains in the singly ionized state, 
C^, up to Fx ~ 5x10^ erg cm“^ s"'. 


= 10 cm^^ and n(H^) = «(e). Fig. 13 a) is a plot of the fractional 
abundances of oxygen ions, 0° to It can be seen that proton 
charge exchange keeps oxygen ionized even in the absence of X- 
rays and that neutral oxygen is essentially absent. As the X-ray 
flux increases increases at the expense of and for Fx > 
fewxlO^ erg cm“^ s ' all the oxygen is in the highest ionization 
state considered here, O"'^. Fig.|5lb) and (c) plot the distribution 
of nitrogen ions, N" to usi ng the charge exchange rate co¬ 
efficients of iLin et all (l2005h and [Kingdon & Ferlandlll99^ . re¬ 
spectively. Their distribution is somewhat different than oxygen 
at low X-ray fluxes because the temperature is not high enough 
for efficient charge exchange by to keep all the nitrogen ion¬ 
ized. Instead the n itrogen is only partially io nized, but more so 
using the results of lKingdon & Ferlandl(ll996l) . As the X-ray flux 
increases it ionizes the nitrogen and increases, and then at 
larger fluxes the distribution starts to resemble those of carbon 
and oxygen. For Fx > fewxlO^ erg cm“^ s * all the nitrogen 
is essentially in the form of (the highest ionization state in¬ 
cluded in our models). 

We expect the fine-structure far-infrared lines from N^, O", 
and to come primarily from the WIM or H ii regions because 
their abundances will be much smaller in neutral clouds because 
of the rapid charge exchange and ion molecule reactions of 
and with H and H 2 . In the WIM, if the dominant processes 
are proton charge exchange. X-ray photoionization, and electron 
radiative and di-electronic recombination, then the distribution 
of lower ionized states of nitrogen and oxygen in the WIM will 
be shifted towards higher ionization states for large X-ray fluxes. 
We expect that the emission from lines such as [O iii] at 52 jum 
and 88 jum, [Nii] at 122 jum and 205 /tm, and [Niii] at 57 /rm, 
will also be suppressed in the presence of a strong X-ray flux, as 
is observed in a number of starburst and active galactic sources. 


for -I- H is fast, ~10 ® cm^ s ', and, while that for -H H is 
-10“'^ cm^ s“' at 8000 K, it is still faster than that for C'^. 


units ~ 224K) that the protons in the warm and hot ISM can ef¬ 
ficiently charge exchange with O, but in the case of nitrogen the 
energy difference is so large (~10,980 K) that the ionized gas 
must be > fewxlO^ K for efficient ionization of N. Furthermore, 
cal culations of the cross sections for -i- N —> -t H dif¬ 

fer dKingdon & Ferlandll996l : lLin et al.l2005h resulting in some¬ 
what different reaction rat es at the temperatures of interest in the 
WIM (lLangeretalJl2015h . 

In ISM gas with atomic hydrogen, the ions will neutralized 
by charge exchange with H because the reaction rate coefficient 

T 3 - ■ 

cm s 

slow. 

Furthermore, in contrast to the charge exchange cross sec¬ 
tions for and with H are fast, ~10“'' cm^ s“' at 8000 K, 
so higher ionization states of N and O will be difficult to sustain 
in gas containing atomic hydrogen. It will be even more diffi¬ 
cult to have an abundance of and in molecular hydrogen 
gas as their reaction rate coefficients with H 2 are very large and, 
indeed, are primary pathways to forming oxygen- and nitrogen¬ 
bearing molecules in shielded regions. Therefore, highly ionized 
states of nitrogen and oxygen are likely to be found only in the 
warm highly ionized gas. 

We have modeled the distribution of the ionized states of ni¬ 
trogen and oxygen under conditions of a warm dense ionized 
medium, similar to the model for ca rbon in Section 13.21 We use 
the photoionization cross-sections in IVerner et al. ^^9^, charg e 
exchange reactio n rat e coefficients from iMcElrov et al.l (1201 31) . 
iLin et alJ (l2005l) . and iKingdon & Ferlandl ( 1996h. and e l ectron 
rec ombination rate coefficients fr om Nahar & PradhanI (Il997h 
and iBadnell et al] (l2003h : [Badnelll (l2006l) . We show in Fig.ISlre- 
sults for nitrogen and oxygen in a dense ionized gas with n(e) 


4. Discussion 

The solutions in Section [3] show that it is possible for a high X- 
ray flux in galactic nuclei to alter the carbon ionization balance 
and reduce the abundance and correspondingly the [C ii] lu¬ 
minosity. In the presence of a high flux of soft X-rays above 
1 keV, a condition encountered in ma ny active galactic nuclei 
(IStacev et al.ll20ir)l : lEbrero et al.ll200^ . the abundance of singly 
ionized carbon is reduced and converted into higher ionization 
states. This reduction occurs primarily in the hot highly ionized 
gas that Alls most of the galactic central zone, to some degree 
in the dense ionized skins surrounding molecular clouds, less so 
in diffuse atomic hydrogen clouds, and very little, if at all, in the 
dense PDRs at the edge of the CO molecular cores. Furthermore, 
in galactic nuclei with a high X-ray flux we would expect t o see 
an inc rease in the dust temperature and infrared luminosity dVoid 
Il99lh . Therefore in galactic nuclei, and in particular in active 
galactic nuclei, we expect a reduction in [C ii] emission relative 
to FIR/IR emission depending on the relative contribution of dif¬ 
ferent ISM components to the [C ii] luminosity. Unfortunately 
the next most abundant ion, does not have fine structure far- 
IR emission lines as its 2s^ ground state has spin zero, and 
although it has spin angular momentum due to its unpaired elec¬ 
tron in the 2S level, does not have a nuclear spin or orbital angu¬ 
lar momentum to break the degeneracy of the two electron spin 
states. Instead to test the effect of X-ray ionization on the carbon 
balance we would need studies of their UV emission. For exam¬ 
ple the [C iv] UV resonance lines are detected in extragalactic 
sourc es and have been used to trace the formation rat e of massive 
stars (iLeitherer & Lamerllll991t [Robert et al1ll993l) . [Civ] UV 
absorption lines have also been used to study the properties of 






















































Langer and Pineda: [C ii] emission from galactic nuclei in the presence of X-rays 



1.0 


0.9 


Nitrogen 


t 

0 

i 

0.8 

- 

n(e) = 10 cm'^ 

T = 8000 K 


1 

! c 

1 

0.7 

- 


1 

- 

0.6 

- 


1 

- 

0.5 

- 


/ 

- 

0.4 

- 


1 

\ 1 
\ 1 
\ 1 

- 

0.3 

- 


X ^ 

- 

0.2 

- f(N) 

-f(N+) 


1 /A « 
v7 \ *' 


1(N2+) 
. f(N3+) 

A ■' 


1 

0.1 

-f(N4+) 

/ \ * 

* \ \ 






\ 


10’“ lO"' 10'' 10° 10’ 10° 10° 10* 

(erg/cm^/s) 


Fig. 5. Fraction of oxygen ionization states (a) and nitrogen ion¬ 
ization states (b and c) versus X-ray luminosity flux, Fx, for a 
dense WIM with n(e) =10 and Ti^=8 000 K. Panel (b) as¬ 
sumes t he charge transfer rates of Lin et alJ (l2005h and panel (c) 
those of iKingdon & Ferlaiidl (1199^ - see text. 


the Galactic Halo (ISavage et alJ200(il) where c arbon is presumed 
to be collisionally ionized in hot coronal gas (iGnat & Sternberg! 

|^o3). 

To illustrate the potential effects of X-rays in reducing [C ii] 
emission we have run a number of simple models of [C ii] lu¬ 
minosity as a function of X-ray luminosity. One model uses the 
ISM conditions of the CMZ and the other the higher densities 
appropriate to AGNs. To model any actual nucleus we would 
need to know many parameters including the properties of var¬ 
ious ISM components - their densities, temperatures and fill¬ 
ing factors, and the distribution and luminosities of the X-ray 
sources. As will be seen below the primary effect of the X-rays is 
to reduce the [C ii] emission from the ionized gas. Thus whether 
X-ray photoionization plays an important role or not depends on 
the relative contribution of the WIM and PDR components in 
galactic nuclei. In the remainder of this section we present sim¬ 
ple models for the contribution to [C ii] emission from the ISM 
components in galactic nuclei and parameterize the effects of X- 
ray ionization on the reduction in [C ii] emission. In one we will 
be guided by the properties of the ISM in the CMZ, but consider 
much larger X-ray luminosities than found there. In the other 
we consider much higher electron densities for the WIM and H 2 
densities in PDRs, as well as a larger volume. Our intent is not 
to model any particular galactic nucleus but to highlight the gen¬ 
eral behavior of the [C ii] emissivity in AGNs and similar active 
nuclei as a function of X-ray luminosity. 


4.1. [Cii] emissivity from X-ray dominated nuciei 

We can estimate the [Cii] luminosity by summing the contri¬ 
butions of the different ISM components over their respective 
volumes. To simplify our model calculations we assume that the 
physical parameters within each of these components are uni¬ 
form and represent typical gas conditions. The total galactic lu¬ 
minosity is 

Lc([C II]) = Vg X n]) erg s-' (18) 

j 


where Vc is the total volume of all [C ii] emitting regions, tjj is 
the fractional volume occupied by each ISM component j, and 
e/([Cii]) is the [Cii] emissivity in erg cm“^ s“' of component j. 
The emissivity for optically thin lines, 

G([Cn]) = (erg cm“^ s“') (19) 

where A 3 / 2 U /2 is the Einstein spontaneous radiation coefficient, 
V 3 / 2 /,i /2 is the transition frequency, and n 3 / 2 j(C^) is the number 
density of the upper ^^ 3/2 state in ISM component j. In the opti¬ 
cally thin limit one can solve exactly for the density of the ^^ 3/2 
state as a fu nction of the density and temperature of the collision 
partner (see [Goldsmith et ani2012l : lLanger et alJl20145) . 


«3/2j(C^) 


2nj(X)n,(C+)e-'^^/^‘“” 
nj.cAX) -F nj{X){l + 2 e-A£/rki„) 


( 20 ) 


where AF =91.25K, n/X) is the density of the dominant col¬ 
lision partner X (= e, H, H 2 ) in each ISM c omponent /, and 
n ;.cr(X) is the corresponding critical de nsity (see lGoldsmith et ^ 
l2012t IWiesenfeld & Goldsmithll2014l) . In this paper we use the 
exact solution, but note that in the limit where the density 
of the collision partner, n{X), is much less than the critical 
density, n„{X), Equation |20] can be simplified to - 

2(n(X)/ncr(X))e-®‘'25^^“”. 
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We define the average emissivity as <e([Cii])> = Lg([Cii])/Vg^ 
which upon substituting the line parameters for [C ii] yields, 

< e([C II]) >= 3.1 X 10-2° ^ njfjP 3 / 2 ,j(C^)xrot,j(C)nj(X) (21) 

J 

where fj is the fractional abundance of with respect to 
all the charge states of ionized carbon for ISM component j, 
P 3 / 2 ,jiC^) - « 3 / 2 j(C^)/nf(C^) is the fractional population of the 
2^3/2 State, XtotjiC) is the total fractional abundance of carbon in 
all charge states with respect to the total hydrogen density (n,- 
n(H) -H 2 «(H 2 ) -I- n(H^)), and nj{X) is the dominant collision 
partner in component j. 

4.2. Galactic nucleus with CMZ model parameters 

In this section we set up a model galactic nucleus to study the 
[C ii] emissivity as a function of X-ray luminosity based on ISM 
parameters of the CMZ. To model the [Cii] emissivity from 
galactic nuclei we need to know the hlling factor, density, tem¬ 
perature, and fractional abundance of for each ISM com¬ 
ponent. We will use the conditions of the ISM in the CMZ as 
a guide for many, but not all components, and will consider a 
much wider range of X-ray luminosity than is found there. The 
dominant ISM components in the nucleus are: 1) a dense warm 
ionized medium (WIM), 2) a low density hot ionized medium 
(HIM), 3) cold diffuse atomic hydrogen clouds (CNM), 4) a 
dense ionized skin (DIS) surrounding molecular clouds, 5) the 
dense PDRs of GMCs, and 6 ) compact H ii regions. We will ne¬ 
glect the contribution from the hot ionized medium because its 
densities are too low, n(e) ~ fewxl 0“2 cm^^ ^nd the compact 
HII regions which, while bright, have too small a hlling factor, 
to contribute to the bulk of the [C ii] emission. We consider the 
galactic nucleus to be a disk with a size similar to the CMZ. We 
assume a radius r - 200 pc which is the boundary of massive 
giant molecular c louds at I ~ 358?7 in the CMZ as traced by CO 
dOka et anil998h and a scale height z = 30 pc (disk thickness 
of 60 pc) slightly lar ger than the scale height of the dense WIM 
dFerri^re et alJ[200^ l4 

We adopt the following densities and temperatures for the 
key ISM components that contribute to the [C ii] emission from 
galactic nuclei: 1 ) the low density warm ionized medium is as¬ 
sumed toJiave_7jdn_=8000Kai^an_^erage density n(e) ~ 10 
cm 2 dCordes & Laziol120031: lRovll201^ : 2) diffuse atomic hy¬ 
drogen clouds with n(H) = 100 cm-2 and Tkin = lOOK, 3) dense 
PDRs with n(H 2 ) = 10° cm-2 and Tkin = lOOK, and for the 
dense ionized skins we adopt n(e) = 20 cm ~2 from the s tudy of 
two clouds at the edge of the CMZ dLanger et al.ll2015l) . These 
values are summarized in Table |4] We also assume a fractional 
abundance of carbon appropriate to the inner Galaxy, x^otjlC) - 
5x10“^, which is an extrapolati on of the metal abundance profile 
in the disk to inside 3 kpc (see lWolfire et al.ll2003l: iPineda et al.l 
120131) . 

We estimate the hlling factor, 77 , of each component from the 
mass and density of the components. The mass of the GMCs 
in the CMZ is estimated to lie in th e range (2 - 6)xl02 Mq 
dOka et al.l Il998t iFerriere et al.l l2007h and we adopt the inter¬ 
mediate value 4x10^ Mq. Their average H 2 densities are in the 
range 10 ^ - 10 ^ cm “2 \\ie estimate the volume hlled by the 
GMCs by dividing the total mass in GMCs, by the density. For 

^ The scale height of the WIM adopted for the CMZ is less than that 
derived from [C iil for the Galactic disk where th e density is much lower 
dLanger et al.ll2014^ : IVelusamv & LangeHl2014h . 


n(H 2 ) ~ 10"^ cm-2, the volume hlled by clouds is ~8xl0‘2 pc^, 
which yields a hlling factor 77 H 2 about 1.1 % of the volume of the 
CMZ, for 10^ cm 2 the volume is only ~8xl02 pc2, or ~ 
0.11% of the CMZ volume. The [Cii] emission from the dense 
PDR layer outside the CO cloud can be estimated from models 
of massive GMCs (IWolhre et al.ll2010|) who hnd that the PDR 
layers lie in the range 1 to 1.5 magnitudes thick regardless of the 
cloud mass, UV radiation held, or hydrogen density. Therefore 
the fractional volume hlled by the PDR is just the difference be¬ 
tween the volume occupied by the cloud out to the edge of the 
PDR minus that of the dense CO core. For example, assuming 
a typical cloud mass of 10^ Mq and a core density of 10^ cm“2 
and PDR density of 102 cm-2, the PDR occupies ~30 percent 
of the total cloud volume, which corresponds to ~3xl0 2 of the 
volume of the CMZ. For more massive clouds, 10° Mq, and in¬ 
terior and PDR densities, 10° cm 2 and lO'* cm 2 , respectively, 
the PDRs occupy only ~7xl0-"^ of the CMZ. 

Here we will consider a range of PDR hlling factors, //(PDR) 
= 1x10-2 to 3x10 2 . The choice of PDR hlling factor has an 
inhuence on the models for two reasons, hrst, as shown above, 
it takes a very large X-ray Hux to alter the carbon ionization 
balance in H 2 gas, and second, if the X-ray Hux is large enough to 
alter the carbon ion abundance it is also large enough to destroy 
H 2 . Thus increasing the X-ray hux only pushes the PDR-like 
layer deeper into the cloud, but does not destroy it. To simulate 
this situation in our model we assume that X-ray photoionization 
does not alter the C^ abundance in the PDR component. 

The dense ionized skins occupy about 30% of the volume of 
the CO core of the GMCs, for a hlling factor //(DIS) ~2xl0 2 , if 
the size of this r egion is similar to th at derived for clouds in the 
region of Sgr E (iLanger et al.ll2015l) . The diffuse atomic clouds 
have an estimated mass ~5xl02 Mq (IFerriere et al.ll200^ which 
corresponds to a hlling factor ~0.03 at a density //(H) = 100 
cm 2 . The ionized gas hlls roughly 95% of the volume of the 
galactic nucleus, of which 80% is t he dense WIM and ~15% 
the hot ionized medium (HIM) (see IFerriere et al.ll2007l) with 
negligible [C ii] emission. 


Table 4. Galactic nucleus model parameters 


ISM 

Component 

n(e) 

cm“^ 

n(H) 

cm“^ 

n(H2) 

cm“^ 

Tkin 

K 

rf 

Mass 

Mo 

HIM 

0.01 

0 

0 

10 " 

0.15 

3x10" 

WIM 

10 

0 

0 

8000 

0.80 

10 ° 

HI clouds 

0.1 

100 

0 

100 

0.03 

5x10" 

CO clouds 

0 

0 

10 '' 

50 

0.0125 

6 x10" 

PDR 

3.5 

0 

10 ^ 

100 

0 .002'’ 

2 x10° 

DIS 

20 

0 

0 

8000 

0.003 

10 '* 


a. The volume hlling factor rj is dimensionless, and the sum for the 
various does not add up exactly to one because of roundoff, b. We 
consider a range of values from 0.001 to 0.003, with 0.002 in the Table 
representing an intermediate value. 


4.2.1. [Cii] luminosity in the central X-ray model 

To demonstrate the effects of X-ray photoionization on the [C ii] 
emission we calculate the [C ii] luminosity as a function of X-ray 
luminosity for a central X-ray source, such as might be generated 
by a massive accreting black hole using the ISM parameters in 
Table |4] We calculate the X-ray Hux as a function of radius as¬ 
suming no attenuation and r 2 . Fig. shows the total [Cii] lu¬ 
minosity, L([Cii]) in units of solar luminosity, Lq, as a function 
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of X-ray luminosity, L(X-ray), ranging from 10^* to 10^^ erg s“' 
(the upper limit is a factor of a few beyond what is observed 
in AGNs). Results are shown for three different choices of fill¬ 
ing factor 77(PDR) = 1,2, and 3x10“^, and it can be seen that 
it takes an X-ray luminosity >10^*^ erg s“' to have a measurable 
effect on the balance and decrease the [C ii] luminosity in the 
galactic nucleus, and >10"^^ erg s“* to have a significant effect. 
At the highest X-ray luminosities it is possible to suppress, by 
an order of magnitude or more, the [C ii] emission from the X- 
ray region of an AGN. The asymptotic value is the contribution 
mainly from the PDRs. 

For a central X-ray source the contributions to the total lu¬ 
minosity will change as a function of galactic radius. In Fig. [T] 
we plot the contribution of the individual ISM components to 
L([Cii]) as a function of distance from the central source, r(pc), 
for two values of the central X-ray luminosity, L(X-ray) = 10"^^ 
and 10"^"^ erg s ', adopting the intermediate filling factor 7;(PDR) 
= 2x10^^. We chose these luminosities as examples because it 
is in this range of L(X-ray) that the photoionization of the WIM 
begins to make a difference in the relative contributions of the 
PDRs and the WIM. For L(X-ray) = 10''^ erg s ' (Fig. [T] top 
panel) the WIM dominates the contribution to [C ii] throughout 
the nucleus, followed by emission from the dense PDRs, while 
the diffuse H i clouds and dense ionized skins (DIS) around the 
GMCs make negligible contributions to L([C ii]). At higher lumi¬ 
nosities the X-rays have a more profound effect on the total emis¬ 
sion, as shown in Fig. [T] (bottom panel), and there is a different 
mix of ISM components contributing to L([C ii]). For L(X-ray) = 
10"'^ erg s ' PDRs dominate the [C ii] luminosity for r < 100 pc, 
while the WIM dominates [Cii] luminosity outside 100 pc, and 
the Hi and DIS components are, again, negligible. Overall, the 
total [C ii] luminosity decreases with increasing L(X-ray) as was 
shown in Fig.|6] The ISM physical conditions in the WIM, DIS, 
dense PDRs, and atomic hydrogen clouds may be different in the 
AGNs, LIRGs, or ULIRGs from those found in the CMZ, but the 
general dependence of the [C ii] luminosity on X-ray luminosity 
will be similar. 


4.2.2. [Cii] luminosity in the uniform X-ray flux model 


Whereas in some galactic nuclei the X-ray luminosity is dom¬ 
inated by a central source, in others the luminosity may result 
from diffuse emission and/or numerous discrete sources spread 
throughout. To model this situation we assume a uniform flux 
throughout the nucleus and calculate the luminosity, L(X-ray) = 
FxA, where A is the area of the boundary of the galactic nu¬ 
cleus. We adopt the same disk model with a radius of 200 pc 
and 60 pc thick used in the CMZ model with a central X-ray 
source. The exact dimensions are not critical because the results 
can be extended to galaxies with X-ray regions of different size 
through the relationship between flux, emissivity, and luminos¬ 
ity for a uniform environment. The ISM parameters are again 
given in Tableland we adopt the intermediate value ? 7 (PDR) = 
2x10“^. The X-ray luminosity for a uniform flux is just propor¬ 
tional to the area of the outermost volume, L(X-ray)=F'xAAGN, 
where Aagn is the area of the AGN’s uniform X-ray source, 
while the [C ii] luminosity is proportional to the emissivity times 
the volume, L([Cii]) = e([C ii])yAGN, where Vagn is the corre¬ 
sponding volume of the uniform X-ray source. Therefore, the 



Fig. 6. The [Cii] luminosity in units of solar luminosity as a 
function of the luminosity of a central X-ray source. The as¬ 
sumed ISM parameters and dimensions of the nucleus are given 
in the text. The three curves represent different values of the fill¬ 
ing factor for PDR and PDR-like conditions in the nucleus. 


results plotted below using the dimensions of the CMZ can be 
converted to different AGN volumes as follows, 

Taga'([Cii]) = (Vagn/V cMz)LcMz{[Cn]) 

( 22 ) 

TACv(X-ray) = (AAcv/^CMz)TcMz(X-ray), 

in which one scales a given AGN X-ray luminosity by the ra¬ 
tio of areas to the CMZ and then determines the CMZ [C ii] lu¬ 
minosity and scales it by the ratio of the volumes to derive the 
AGN [C ii] luminosity. In Fig. [8] we plot the [C ii] luminosity 
arising from a galactic nucleus with uniform X-ray source. The 
luminosities of the individual ISM components are plotted as a 
function of X-ray luminosity. The dense warm ionized medium 
dominates the [Cii] luminosity up to L(X-ray) ~10"'"' erg s“', 
after which PDRs, with a filling factor /^(PDR) = 2x10“^, domi¬ 
nate the emission. Finally, in Fig.|9]we compare the [Cii] lumi¬ 
nosity of the central and uniform X-ray models as a function of 
L(X-ray) assuming the same CMZ ISM conditions and n(e)=10 
cm“^ for the WIM and the same value for //(PDR) = 2x10“^. The 
two solutions for L([Cii]) are very similar. At low X-ray lumi¬ 
nosity the flux is too low to influence the C^ abundance and the 
emissivities are independent of L(X-ray), while high luminosi¬ 
ties the C^ abundances are reduced in the WIM throughout the 
volume and the PDRs dominate the emission. At intermediate X- 
ray luminosities the solutions for the central source reduces the 
C^ abundance primarily in the interior (smaller r) while having 
less of an effect near the outer part of the nucleus, and this tends 
to average out to a solution close to that of having a uniform flux. 
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Fig. 7. The [C ii] luminosities of the individual ISM components 
are plotted as a function of radius in a galactic nucleus with a 
central X-ray source with L(X-ray) = 10"^^ (top panel) and 10“*"^ 
erg s ' (bottom panel) with a filling factor //(PDR) = 2xl0~^. 
The assumed ISM parameters and dimensions of the nucleus are 
given in the text. The dense WIM dominates the [C ii] emission 
at L(X-ray) =10""' erg s'* but is increasingly suppressed with 
increasing L(X-ray) until the PDRs dominate the emission. 

4.3. AGN representative model of [Cii] emission 

In the CMZ model we adopted an average electron density n(e) 
= 10 cm'^, but the most central region has been est imated to 
have a much higher density (c .f. iFerriere et alJl2()()^ . For ex¬ 
ample, iMezger & Pauls! (Il979ll used thermal radio continuum 
radio observations to derive a central n(e) ~ 26 cm'^ while 
iMehringer et al.l (119921 Il993h derive even higher electron den¬ 
sities outside compact HII regions near Sgr B1 and B2 and de¬ 
rive an average electron density of ~ 60 and 80 cm'^ over two 
regions ~ 700 and 170 pc^, respectively. The conditions encoun¬ 
tered in the most energetic region of the CMZ is likely to prevail 
over a much larger volume in AGNs. In addition, the [N iii] lumi- 



Fig. 8 . The luminosities of the individual ISM components aris¬ 
ing from a galactic nucleus with a uniform X-ray source, and 
physical parameters given in the text and Table 0] are plotted as 
a function of X-ray luminosity. 



Log (L{X-ray)) (erg/s) 


Fig. 9. The [C ii] luminosity for the central X-ray source (black 
line) is compared to that of the uniform flux model (red line). 
The ISM properties are the same in both models. 


nosity in some galactic nuclei is comparable to that of [N ii] but 
requires higher electron densities to excite its ^^ 3/2 level. Thus 
it is reasonable to consider a model for AGNs with much higher 
WIM densities tha n the CMZ. In add ition, the PDR densities are 
also likely higher (IStacev et al.ll2010h . To illustrate the impact of 
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Log {L(X-ray)) (erg/s) 


Fig. 10. The fraction of singly ionized carbon in the WIM versus 
X-ray luminosity in the uniform X-ray flux model, for a range of 
WIM electron densities. 


higher electron densities we first demonstrate the sensitivity of 
the WIM [C It] luminosity as a function of n(e) on a stand alone 
basis and then in a simple model of AGNs. 


4.3.1. [Cii] luminosity from the WIM 

In the model adopted here for the [C ii] luminosity from active 
galactic nuclei the WIM dominates the emission of [C ii] over a 
wide range of X-ray luminosities. However, the model adopts a 
fixed electron density in the WIM, n(e) = 10 cm“^. To under¬ 
stand better the sensitivity to the choice of this parameter we 
have also calculated the luminosity from the WIM, Lwim([C ii]), 
versus the electron density. Fig. [10] shows the fraction of singly 
ionized carbon in the WIM as a function of X-ray luminosity 
for a range of WIM electron abundances, n(e) = 1 to 25 cm“^. 
As expected, the smaller the electron abundance the more pro¬ 
found is the effect of X-rays on /(C^) for a given X-ray lumi¬ 
nosity, because the electron recombination rate of cannot as 
efficiently offset the X-ray photoionization rate. In the absence 
of X-rays the luminosity from the WIM would be expected to 
decrease roughly like n(e)“^ below the electron critical density 
{ncr{s)~ 40 cm“^). However, we expect that the [C ii] luminosity 
from the WIM decreases faster than the density n(e) in the pres¬ 
ence of a large X-ray flux because not only does the collisional 
excitation rate decrease, but so does the fraction of available C^. 
This nonlinear effect can be seen in Fig. [TT] which plots Lwim 
as a function of L(X-ray) for the range of n(e) plotted in Fig.fTOl 
In Fig. [T2| we plot the total [C ii] luminosity, Ltot as a function 
of L(X-ray) for a range of n(e) in the WIM, and it can be seen 
that the drop in [C ii] luminosity can be quite large if the WIM 
electron density is higher, such as might be expected in active 
galactic nuclei. 



Fig. 11. The [C ii] luminosity arising from the WIM versus X- 
ray luminosity, for a range of electron densities in the WIM. The 
remaining ISM parameters are the same as that for the uniform 
model. 



Log {L{X-ray)) (erg/s) 


Fig. 12. The total [C ii] luminosity (all ISM components) versus 
X-ray luminosity, for a range of electron densities in the WIM. 
The remaining ISM parameters are the same as that for the uni¬ 
form model. 


4.3.2. [Cii] luminosity from an AGN 

We model a hypothetical AGN assuming a disk geometry with 
a central X-ray source. The WIM electron density is n(e) = 50 
cm“^ and the PDR has n(H 2 ) = 5x10^ cm“^. The molecular mass 
in active galactic nuclei is larger than the CMZ but as the density 
of the giant molecular clouds is an order of magnitude higher the 
filling factor is not much different from what we calculated for 
the CMZ and we assume 7;=0.002. One other difference is the 
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Fig. 13. The total [C ii] luminosity of a hypothetical AGN versus 
galactic radius for a range of X-ray luminosities. The key labels 
the log of the X-ray luminosity and ranges from LogL(X-ray) 
= 40 to 48. The decrease in [Cii] luminosity with increasing 
L(X-ray is due to the destruction of in the WIM close to the 
source of X-rays. The radius over which the X-rays photoionize 
increases with increasing X-ray luminosity. 


scale height of the gas which is likely to be higher than in the 
CMZ due to a higher pressure and we adopt a scale height for 
CO clouds and the WIM of 100 pc and 200 pc, respectively. The 
[C ii] luminosity for this AGN model is plotted in Fig. [13] as a 
function of galactic radius from the central source for a range 
of L(X-ray). The X-ray flux close to the central accretion hole 
ionizes the and reduces the luminosity of the WIM as can be 
seen in the decrease in the luminosity as a function of X-ray lu¬ 
minosity. As one goes further from the central source the WIM 
increases and contributes to the total. At a few hundred pc there 
is roughly an order of magnitude difference in the total luminos¬ 
ity. The T([C ii]) dependence on L(X-ray) is similar to the CMZ 
example except that the luminosity is larger due to higher elec¬ 
tron and H 2 densities. 


4.4. Observational test case 


In Fig.[T4| we compare our model calculations to a sample of six 
extragalactic sources for whic h the [Cii], FIR, a nd X-ray lumi¬ 
nosities have been measured dStacev et al.ll201(]l and references 
therein). We al s o incl ude two high-redshift AGNs compiled by 
iGullberg et aki (l20L5h in [C ii] and FIR and by IVignali et al.l 
( 2005h and lLaw et al.l (l2004h in X-rays. We calculated the model 
[C ii]/FIR ratio by dividing the predicted [C ii] luminosity for dif¬ 
ferent PDR filling factors (Fig. |6ll by the FIR lu minosity of the 
CMZ, Lfir = 4 X 10* Lq (iLaunhardt et al.ll2002h . The observed 
X-ray luminosities cover 2 to 10 keV and, because the spectral 
energy distribution increases at lower energies the actual lumi- 



Fig. 14. The predicted [C ii]/FIR luminosity ratio versus X-ray 
luminosity for different PDR filling factors (see Fig. |6ll. Also 
plotted are data points for six AGNs (red b ullets) observed 
in [Cii], FIR, and X-rays (IStacev et al.ll2010l) and two high- 
re dshift AGNs (black sq uares) ob served in [C i i ] (as c ompiled 
in iGullbe rg et alJ ( 20151)') . FIR by IVignali et akl OOOSh . and X- 
ravs bv iLawetal. (12004). We assumed 30% uncertainties for 
the [C ii]/FIR ratios from Stacev et"^ d2010h and for the X-ray 
luminosities. 


nosity would be a factor ~3.8 larger for E> 1 keV assuming the 
spectral energy distribution in Equation [TSl with F = 1.9. Thus 
the data points would shift slightly to the right by ~0.59 on the 
log scale, which is still consistent with our model calculation. 
We see that for this sample of AGNs the [C ii]/FIR luminosity 
ratio is in good agreement with our model predictions, showing 
significantly lower [C ii]/FIR ratios for L(X-ray)> 10"^* erg s“'. 

4.5. Central Molecular Zone [Cii] luminosity 

The Galactic Central Molecular Zone provides the best opportu¬ 
nity to check our model of [C ii] luminosity in the presence of 
X-rays because we have the best information about the state of 
its ISM components, at least as compared to external galaxies. 
However, in practice the uncertainties in its ISM conditions are 
too large to allow more than an estimate of the [C ii] luminosity. 

Unfortunately, the Galactic Central Molecular Zone is not a 
strong source of X-ray luminosity. Chandra has det ected >9000 
discre te stellar sources in the innermost CMZ (iMuno et al.1 
I 2 OO 9 I) . but most of these have luminosities in the range L(X- 
ray) = 10*' and 10** erg s“* with a few having maximum values 
~ fewxlO*"' erg s“*. Adopting an average luminosity per source 
of 10** erg s“' yields a maximum extended luminosity of order 
10*® erg s ', which is too low to have much of an effect on the 
[C ii] luminosity based on the results in Fig. [8] The strongest 
point source in the CMZ, IE 1740.7-2942, is believed to be 
an accreting stellar bl ack hole and has a lu minosity, L(X-ray ~ 
fewxlO** erg s ' (see lMalonev et al.lll997h . The black hole as¬ 
sumed associated with S^r A at the Galactic Cen ter has a current 
luminosity ~10** erg s“' (iBaganoff et al.1120011) . but may have 
been brighter in the recent past. Thus X-rays likely have little 
effect on the [Cii] luminosity of the CMZ. and so it is not a 
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good source for testing the effects of X-rays on [C ii] luminosity. 
Nonetheless, the CMZ [C ii] luminosity does provide an oppor¬ 
tunity to validate whether our model predicts the right order of 
magnitude for the [C ii] luminosity. 

To calculate the [C ii] luminosity for the CMZ we use the 
CMZ physical dimensions adopted in Section 14.21 and ISM pa¬ 
rameters in Table |4]for the atomic clouds, PDRs, and DIS com¬ 
ponents, butmodify_AeWIM model to ht with the WIM electron 
prohle (iFerriere et al.l 1200^ . as discussed below. The atomic, 
molecular, and ionized gas in the CMZ is concentrated at the 
center and decreases with radius and distance above and be¬ 
low the plane. We are interested in calculating the [C ii] lumi¬ 
nosity arising from the volume of the CMZ so the distribution 
of the atomic clouds and the PDRs and DISs associated with 
the giant molecular clouds is not important under the assump¬ 
tion that the clouds have the same characteristics throughout 
the CMZ volume, as assumed here, and as long as the X-ray 
flux has a minimal impact on the C^ abundance. The WIM den¬ 
sity, in contrast, does change throug hout the CMZ. IFerriere et al.l 
( 2()()7|) summari z es the result s from iLazio & CordesI (Il998l) and 
ICordes & Lazio! (l2002l l2003h and expresses the electron space 
averaged density as the sum of three terms, one of which, < 
n(e )>3 (r, z), dominates in the CMZ with a prohle that decreases 
as a Gaussian in r and z (the other two terms are smaller in ab¬ 
solute value but have much larger scale factors and so decrease 
very slowly in the CMZ). We have simplihed this expression by 
neglecting the offset terms in < n(e)> and setting < n(e) >i and 

< n(e )>2 constant in the CMZ, which yields 

< n(e) > (r,z) - -H 0.14 cm“^ (23) 

where the peak density is ~ 10 cm“^ at the center and decreases 
with Gaussian parameters ro = 145 pc and zo - 26 pc. 

We have calculated the [C ii] emissivity throughout the CMZ 
disk (r < 200 pc and |z| < 30 pc) using the model of the CMZ 
discussed above and assuming a central source with X-ray lumi¬ 
nosity 10^^ erg s“'. We integrate the emissivity over the volume 
of the CMZ to derive an emergent [C ii] luminosity, L([C ii]) ~1 
- 9 xlO^ Lq, for ? 7 (PDR) = (1 - 3)xl0“^. This luminosity hardly 
differs from that without an X-ray source because, as shown ear¬ 
lier, the X-ray luminosity is too low to modify the C^ distribu¬ 
tion in the ISM components. For comparison we have calculated 
the [C It] luminosity measured by BICE usin g the line intensity 
contours in Fig. 8 of iNakagawa et akl (Il998h integrated over the 
area of the CMZ. We estimate that BICE observes L([Cii]) ~ 
4.5x10^ erg s“'. Our L([Cii]) model of the CMZ and the BICE 
observations agree within a factor of 2 or better which is proba¬ 
bly about the best that we can expect given the uncertainties in 
ISM parameters assumed in our model calculation. In summary, 
the reasonable agreement of our [C ii] emissivity model and the 
observed [C ii] luminosity of the CMZ lends conhdence to the 
approach presented here. 

5. Summary 

We have explored the influence of X-rays on the ionization bal¬ 
ance of carbon in the interstellar medium with a focus on its 
effect in galactic nuclei with large X-ray luminosity. In the pres¬ 
ence of a large X-ray flux we hnd that C^ can be converted to 
higher ionization states, C^^ with J > 2. The effect of X-rays 
on the [C ii] luminosity of galactic nuclei depends on the X-ray 
luminosity and the details of the properties of the ISM, namely 
the density, temperature, and hlling factors for each of the ISM 
components. The biggest impact is on the ionization distribution 


in the WIM, where C^ is more readily converted to C^^, C^^, 
etc., thus reducing its contribution to the [Cii] luminosity. 

Thus X-ray photoionization is another mechanism to con¬ 
sider in understanding the reduction of the [Cii] to far-IR lu¬ 
minosity ratio in active galactic nuclei. It appears that, on aver¬ 
age, it requires a soft X-ray luminosity greater than about 10"^^ 
- 10"^^ erg s“* to have a noticeable effect on the [Cii] luminos¬ 
ity. Our Galaxy does not fall into this category, but many others 
do, in particular some Luminous Infrared Galaxies (LIRGs) and 
Ultra Luminous Infrared Galaxies (ULIRGs), and many Active 
Galactic Nuclei (AGNs). If the primary X-ray source is an ac¬ 
creting black hole then the [C ii] will be suppressed in its vicin¬ 
ity but not necessarily throughout the entire galactic nucleus. 
However, as discussed above, AGNs have, on average, consid¬ 
erably larger X-ray luminosities with L(X-ray) ranging up to a 
few xlO"^’ erg s^^ (lEbrero et al.ll2009l) . In the AGNs with L(X- 
ray) > 10"^^ erg s ', depending on the details of the ISM param¬ 
eters, X-rays can suppress the [Cii] luminosity. It is precisely 
these sources that consistently have a much lower [C ii] luminos¬ 
ity comparedjojheirfer-infrared dust luminosity (see summary 
bv iGullberg et al.l[2015l) which may be due, in part, to the influ¬ 
ence of X-rays on the photoionization into higher ionization 
states. 

We have also briefly considered the ionization balance of ni¬ 
trogen and oxygen in the dense warm ionized medium under 
the influence of X-rays. We find that the lower ionization states 
of nitrogen and oxygen are suppressed at high X-ray fluxes and 
at the highest flux values these elements are essentially in the 
highest ionization state. We did not calculate detailed models 
of the emission from the far-infrared lines of the lower ioniza¬ 
tion states, but it is apparent, by comparison with the results for 
carbon, then the emission from these tracers will be suppressed 
at high fluxes. [N iii] is observed to be strong in some galaxies 
which requires an ionization source for N and N^. The usual 
explanation invokes an EUV flux, but X-rays are an alternative 
explanation and have the advantage of being present throughout 
galactic nuclei. Thus X-ray photoionization may explain some of 
the decrease in emission from carbon, nitrogen, and oxygen ions 
and the increase in emission of higher ionization states, such as 
[N III] and [Oiii], in active galactic nuclei. 

We have been guided by conditions in the Galactic Central 
Molecular Zone, but AGNs will have different ISM parame¬ 
ters of density, temperature, and hlling factor for the various 
gas components and each one would need to be analyzed on 
a case by case basis. Eor comparison we adopted a model of 
an AGN with higher WIM and PDR densities. We also com¬ 
pared our [C ii] luminosity models with a small data set of extr- 
galactic sources with high X-ray luminosities. We hnd that their 
[C ii] to EIR ratio decreases signihcantly at an X-ray luminosity 
consistent with our model calculations. Einally, to conhrm that 
the [C ii] emissivity model we use is reasonable we calculated 
the [C ii] luminosity for the CMZ where conditions are reason¬ 
ably well known, at least in comparison to extragalactic sources. 
The Galactic CMZ model produces a [Cii] luminosity ~ 1 - 
9 xlO^ Lq. depending on the details of parameters chosen for 
the WIM and PDRs. This [C ii] luminosity compares reasonably 
well, given all the uncertainties in the ISM conditions, to the 
[Cii] luminosity of ~ 4.5x10^ L© we estimate from the BICE 
observations. 

In conclusion. X-ray photoionzation can alter the ionization 
balance of the ionization states of carbon, nitrogen, and oxygen 
primarily in the ionized ISM, resulting in a signihcant fraction of 
their abundance shifting tohighly ionized states in active galactic 
nuclei. Therefore the effects of X-rays on the ionization balance 
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needs to be included in any modeling and interpretation of car¬ 
bon, nitrogen, and oxygen far-IR line emission in active galactic 
nuclei. 
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